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Abstract 

Compatibility conditions between the (global) spectrum of an n-mode Gaussian state and the 
spectra of the individual modes are presented, making optimal use of beam splitter and (two-mode) 
squeezing transformations. An unexpected bye-product of our elementary approach is the result 
that every two-mode Gaussian state is uniquely determined, modulo local transformations, by its 
global spectrum and local spectra - a property shared not even by a pair of qubits. 
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The quantum marerinal problem has attracted considerable interest in quantum informa- 
tion theory Given a multipartite system, it asks: what kind of 
spectra for the subsystem density operators are consistent with a given spectrum for the 
density operator of the full system? The Gaussian quantum marginal problem (detailed be- 



low) has 
earlier 



3een solved recently [10|, lUl] (As noted in Ref . [ll|] , the three-mode case was known 
12|). Our approach to this problem makes effective use of beam splitter and two- 
mode squeezing transformations. In the case of two modes it is shown that every Gaussian 
state is uniquely determined, modulo local canonical transformations, by its global spectrum 
and local spectra; in particular, the entanglement is fully determined by these spectra. 

Consider a Gaussian state of a system of n-modes, represented by density operator p. 
The mean values of the position and momentum variables qj, pj have no role to play in our 
considerations, and so we assume that these mean values vanish. Such a zero-mean Gaussian 
state is fully described by its 2n x 2n covariance matrix V. 

The reduced state pj of the j^^ mode, obtained by tracing out from p all other modes, 
is also a zero-mean Gaussian state. With the phase space variables assumed arranged 
in the order qi,pi; q2,P2', ■ ■ ■ qn,Pn the j^^ 2x2 block along the leading diagonal of V 
represents precisely the covariance matrix of the reduced state pj. Through (independent) 
local canonical transformations G 5*^(2, R) on each mode we make all the 2x2 blocks along 
the diagonal of V multiples of identity. The covariance matrix of the j^^ mode will then be 
of the form diag(mj, mj). It corresponds to a thermal state, with temperature T{mj) which 
is a monotone increasing function of nij. Being thermal, pj has the spectral resolution pj = 
[1 — C,{mj)] Y^'kLo^{^jT^''\^jk) {njk\ ■ The parameter Orrij) is another monotone increasing 
function of m^, and |njfc)'s are the energy eigenstates of the j*^ oscillator. Clearly, the 
eigenvalue spectra of the pj's are determined by, and determine, the local spectral parameters 



Using an appropriate (nonlocal) canonical transformation S G Sp{2n, R) the covariance 
matrix V can be decoupled and brought into the canonical form V^'^^ of independent oscilla- 



tors in thermal states 



131 ]: V^^^ = SVS^ =diag(Ki, ki ; ^2, ^2 ; ■ ■ ■ ; Hn, k^)- The associated 



density operator p^^^ thus has the spectral decomposition 

n cxD 

= m^-a^^)]T.a^^r>,,)M (1) 

j=l fc=0 

Since p^^^ and the original p are unitarily related, the spectrum of p is the same as that 



of p^^\ It is clear that this global spectrum and the n-tuple of global spectral parameters 
K2, - ■ ■ , Kn) determine each other. 

We may now ask what are the constraints connecting the global spectrum of a Gaussian 
state to its local spectra. In view of the invertible relationships just noted this Gaussian 
quantum marginal problem is equivalent to seeking the compatibility constraints between 
the global spectral parameters { kj } and the local spectral parameters { mj }. Interestingly, 
the answer can be given in the form of necessary and sufficient conditions. 

Theorem: Let m = (mi, m2, m^, ■ ■ ■ , m„) and k, = K2, ■ ' ' 5 i^n) be the local and global 
spectral parameters of an n-mode Gaussian state, written in nondecreasing order. These 
are compatible iff 

k k 

/c = 1,2, ■ ■ ■ ,n, (2) 

i=i i=i 

n— 1 n—1 

mn-^rrij < Kn-Y^ nj. (3) 

j=i j=i 

Remarks: What this claim means can be clarified by stating it in two parts. Suppose a 
Gaussian state is given. Its local spectral parameters mi, m2, ■ ■ ■ , m„, and global spectral 
parameters ki, K2, ■ ■ ■ , Kn are certain to meet these inequalities (with /ti > 1). Conversely, 
given a set of local and global spectral parameters meeting these inequalities (with ki > 1), 
we can certainly construct a physical Gaussian state with these parameters. 



The first part of the theorem was essentially proved by Hiroshima [10(] . But the full 
theorem in this form was formulated by Eisert et al. [UL] who presented an inductive proof 
for the second part. Our proof of both parts will be seen to be constructive, consistent 
with the elementary nature of the theorem, and it rests in an essential manner on a fuller 
appreciation of the two-mode situation. 

Given two vectors m, k, & R", we will say k, dominates m if m and k, after their 
components are rearranged in the nondecreasing order, obey the set of + 1 inequali- 
ties (2), (3). This definition is such that permutation of the components of m or k, does 
not affect dominance. Thus (9, 7, 8, 6, 12, 11, 10) is dominated by (5, 2, 18, 4, 1, 12, 3), since 
(1, 2, 3, 4, 5, 12, 18) manifestly dominates (6, 7, 8, 9, 10, 11, 12). Further, dominance so defined 
is transitive: k dominates m, and m dominates m', together imply k, dominates m'. 



In the Schur-Horn case 



14| wherein m corresponds to the diagonal entries of a hermitian 



matrix and k to its eigenvalues, the last inequality in (2) becomes an equality. It is clear 



that (3) is subsumed by (2) in that case. 

The case of two modes: This case is of interest in its own right. Further, it possesses an 
aspect which seems to be unique, not shared by any other system. Finally, our analysis of 
the n-mode case relies critically on repeated applications of the two-mode result. Hence we 
begin with a direct proof of the theorem in the two-mode case. 

Lemma: The parameters mi < m2 and 1 < < ^2 are compatible for two-mode 
Gaussian states iff 



mi m2 > Ki + K2, 
m2 — mi < ^2 — Ki. 



(4) 



Note that the condition mi > ki is subsumed by (4). 

Proof of Lemma: The covariance matrix can, through local unitary (canonical) transfor- 
mation G 5*^(2, R) X S'p(2, i?), be brought to the form 



V 



(5) 



'^mi ^ 

mi kp 
kx m2 

kp 7712 

The global spectral parameters ni, K2 are related to the local mi, m2 through the symplectic 
invariants llSll 



V 



-tr {QVQ^V) = kI + Kl = ml + ml + 2kxk 
detV = {K1K2Y = (mim2 — A;^)(mim2 — /Cp). 



(6) 



These immediately imply 



K1K2 < mim2, 
1^1 + ^1 > m^ + mg, if k^kp > 0, 
k\ + k\ < m\ + m2, if k^kp < 0, 

equality in the first inequality holding if fca, = = kp. These inequalities imply 

1^2 — 1^1 > 1^2 — f^i, when kxkp > 0, 
K2 + Hi < m2 + mi, when k^kp < 0. 



(7) 



(8) 



This much is immediate from the symplectic invariants. What remain to be proved are : 
i<'2 — i^'i > 1^2 — TTii when kxkp < and K2 + k,i> 777.2 + ttii when kxkp > 0. 

To prove these we reinterpret (6) as simultaneous expressions for k^, kp in terms of 

Ki, K2; mi, 7712. 

kxkp = [ + 4) - {ml + ml) ]/2, (9) 
kl + kl = -^—[mlml-Kl4 + klkl]. (10) 

^ ^ 77717772 ^ ^ ^ ^ ^ PJ V / 

It is clear that real solutions for k^ and kp will exist iff ' /c^ + /c^ ' > ' 2| k^kp \ '. That is, iff 

mim2 — \kj.kp\ > K,iK,2. (11) 
With use of (9) for k^kp, this last condition reads 

K2 — Ki > m2 — mi, when k^kp < 0, 

K2 + K,i < 7?72 + 7?7i, whcn kxkp>0. (12) 

Proof of the Lemma is thus complete. 

Two types of simple transformations on any pair of modes characterised by annihilation 
operators aj, deserve particular mention; they play a key role in our proof of the theorem. 
The first, Sq, corresponds to the compact transformations % — > cosOaj + sin^afe, — > 
— sin^Oj + cos^Ofe, and therefore is represented by Sg — cos^ctq (8) ctq + sm9ia2 <S> ctq & 
Sp{4:, R), < 9 < 2%, where uo is the 2x2 unit matrix and a2 is the antisymmetric Pauli 
matrix. Physically, Se is a beam splitter with transmitivity cos^ 6. The second one, S*^, is 
noncompact and corresponds to squeezing transformations aj — > cosh/xa^ + sinh/xa^, — > 
cosh //Ofc+sinh a], and is represented by 5'^ = cosh // cro(H)(7o+sinh (7i(8)(73 e 'S'p(4, i?), < 
II < oo. 

It is easily verified that when the covariance matrix V, Eq. (5), has kp — k^ = k, it can 
be diagonahzed by the beam sphtter transformation V SeVSj' , with 9 fixed through 
tan 2^ = 2k/ {m2 — mi). And k,2 + Hi will precisely equal m2 + mi in this case. Similarly, 
if kp = —kx = k > 0, then V is diagonahzed by the squeezing transformation V —>■ S^V S'^ , 
with tanh 2/1 = 2k/ (7712 + ^1), and one will find K2 — Ki — m2 — mi in this case. 

Conversely, suppose we start with the canonical form V^^^ — diag {ki, ki;k2, K2), and we 
wish to achieve through symplectic congruence V^^^ — > SV^^^S'^, S e 5*^(4, R), a covariance 
matrix with diagonals 7?7i, 7772- If mi < 7772 are such that 7772 < k,2 and K2 + Ki — 777,2 + mi. 



such a redistribution of ki, K2 among mi, m2 can always be achieved through a beam sphtter 
transformation S0. Under Se we have m2 + mi = + and m2 — mi = cos 26 (/«2 — ni)- 
On the other hand, if m2 > k,2 and K2 — Ki = m2 — mi, so that Ki and fi;2 are enhanced 
by equal amounts to mi, m2, this can be achieved through a squeezing transformation 5*^. 
Under S*^ we have m2 — mi = K2 — and m2 + mi = cosh 2/i (k2 + 

Our Lemma is similar to Lemma 5 of Ref. [Ill], but our proof is direct and constructive. 
There is an important distinction in content as well: while theirs claims that m2 — mi = 
K,2 ~ iff ^2 = ^^2 and mi = Ki, we have just demonstrated that if m2 — mi = /t2 — Ki 
then m2 + mi could equal cosh 2/i (^2 + Hi) for any < /i < 00, not just /i = 0. Indeed, 
this distinction is central to Stage 2 of our proof of the second part of the main theorem, 
the part which distinguishes the present symplectic situation from the Schur-Horn case. 

Returning to Eq. (10), if we are given values for the expressions ' + 6^ ' and ' ab ' with 

+ b"^ > 2\ ab\, the solution for (a, 6) is unique [{a,b) and {b,a) are not distinct solutions 
for our purpose]. This innocent looking observation leads to a surprising conclusion. 

Proposition: Specification of the local and global spectra of a two-mode Gaussian state 
determines uniquely the state itself, modulo local unitary transformations. 

States of a pair of qubits share a similarity with two-mode Gaussian states in important 
respects. For instance, positivity under partial transpose is a necessary and sufficient con- 
dition for separability and nondistillability in both cases. But a statement analogous to the 
above proposition is not true for a pair of qubits! 

Proof of main theorem: Assume we are given a (zero-mean) Gaussian state, or equiva- 
lently, an acceptable covariance matrix V, the 2x2 blocks along the leading diagonal of V be- 
ing of the form diag(mj, mj). The global spectral parameters { } are immediately defined 



by V 



13| . It is assumed that m = {mi,m2, ■ ■ ■ mn) and k = [ni, ^2, ■ ■ ■ , «„) are arranged 



in nondecreasing order. Let P^ denote the product k,iK,2 ' ' ' i^n and let P^ = mim2 ■ ■ ■ mn- 
Clearly, = det V < Pm, equality holding iff V is diagonal, i.e., iff mj = kj, j = 1,2, ■ ■ ■ ,n. 
Our task is to prove that k, dominates m. 

Choose a pair 1 < j < k < n such that the 2x2 block (in the off-diagonal location) 
connecting the j^^ and k^^ modes is nonzero. We can arrange (through local rotations) 
this block to be diagonal. Let us 'diagonalize' this 4x4 part of the covariance matrix 
using an appropriate two- mode canonical transformation G Sp{A, R), so that mj and m^ are 
transformed to rhj and rhk respectively, the other diagonal parameters remaining unaffected. 
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It is be noted that the new m dominates the original m. That this is so follows, in the case 
k < n, from the facts rhj < irtj and rfij + fhk < mj + rrik- In the case k = n it follows from 
the additional fact that if is less that it is so by a magnitude which does not exceed 
the magnitude by which rhj is less than rrij {rhk — rhj > rrik — rrij). Further, rhjfhk < rrijmk- 

Denote by m' the new diagonal m-parameters arranged in nondecreasing order by corre- 
spondingly permuting the oscillators. Since rhjfhk < rrijmk we have Pm' < Pm- 

For purpose of clarity, let us carry out this process one more time. The parameters ml 
will then go to m" dominating m', with P^" < Pm'- It follows from the transmitivity of 
dominance that m" dominates m. 

It is now clear that when this process is iterated, m goes through a sequence of interme- 
diate values, the value at every stage dominating the previous value, and correspondingly 
Pm steadily decreasing, until Pm reaches or, equivalently, until V becomes diagonal. This 
completes proof of the first part of the theorem. 



10|. Pm played 



The elementary nature of our proof may be compared with that of Ref. 
the role of 'profit function' monitoring progress of this diagonalization process. 

To prove the second part assume, conversely, that we are given the global and local 
spectral parameters k, m E P". Assume that these are compatible: i.e., k dominates m, 
with Ki > 1. Our task is to construct a Gaussian state with these properties. In other words 
we have to present a canonical transformation 5* G Sp{2n, R) which acting on a covariance 
matrix V = diag(Ki, /ti ; k,2, 1^2] ■ ■ ■ ; i^n, i^n) will produce a covariance matrix SVS'^ with 
the target diagonal values m. We build such an S' as a product of n — 1 specific two-mode 
transformations, evolving m*^°^ = k successively through a sequence of intermediates m^^\ 
rrS^\ ■ ■ ■ to finally m*-"""'^^ = m. It will be manifestly clear that m'-'^-* dominates m^'^'^^^ at 
each stage. For clarity, this process is implemented through four elementary stages. 

Stage 1: Since m*^°) = k dominates m, we have mi > rrif^ = ki. Suppose mi = mi'^ + ei, 
ei > (one will move to the next step if mi = mi'^). Let ji be the least integer 
< n such that m^^^ > mi. Carry out a beam splitter transformation Sg between the 
first and ji^^ mode so that the corresponding diagonal elements {mf\m^^^) get redis- 
tributed to {mi^^ + ei,mj^^ — ei) = {mi,m^^^ — ei), with no change in the other diago- 
nal entries: m*^°) = {mi'\m2'\ ■ ■ ■ , m^'^) m^^^ = {mi, mg^"*, ■ ■ ■ , mj^^ — ei, ■ ■ ■ , m^^^) = 

We can repeat this process. Let m2 = m^2^ + €2- By hypothesis €2 > (this is so even if 
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ji had equalled 2). Assume 62 > (if €2 — 0, one moves to the next step). Let j2 be the 



smallest integer < n such that mfj > m2 [Clearly, j2 can be as small as ji, but not any 
smaller]. Carry out a beam sphtter transformation on the 2'^'* and jf" modes so that the 
corresponding diagonal elements {m^2\''^fj) S^t redistributed to {m2,mj^J — 62) to produce 
m^^), leaving the other diagonals unaffected. 

If we are able to repeat this process only i times we have, at the end of it, 

m(^) = ( mi, ma, ■ ■ ■ , nif, mf}^, mfl^, ■ ■ ■ , m^^^ ), (13) 

with mf^ < ruj, V£+l < j < n — 1, and m^'^ = m!^^ — Kn- What we have done so far 
is identical to what one would have done in the Schur-Horn situation. Clearly, the beam 
splitter transformations carried out so far affected neither the sum of the diagonal entries of 
m^'^ nor its n*^ entry. Consequently, the difference m^*^^ — YJjZl 'nT'f^ has remained the same 
for all < A; < £. 

Stage 2: Define = E]=i mj - E"=i mf. It is clear that 6^''^ = 5^, for A; = 1, 2, • • • , /. 
In the Schur-Horn situation S^^^ vanishes by hypothesis. We will now employ two-mode 
squeezing transformations Sf^ to rectify this 'departure' from the Schur-Horn situation. 

We know that 5^^^ = 5^°^ is nonnegative. Assume 5^°^ > ( if 5^°^ = 0, one will move 
directly to Stage 4, as will become evident below). Define e^+i = m^+i — mfli. Assume 
5'^^^ > 2ei+i ( if this is not the case one will move to Stage 3). Carry out a two-mode 
squeezing transformation S*^ between the (£ + 1)*'^ and n''^ modes, raising the corresponding 
diagonal entries mfli, m^^ — m^^ — Kn by equal magnitude to m^+i, m^^ + e^+i with no 
change in the other diagonal entries, so that 





= (mi, 


(£+1) (e+i)\ 






(^) 
= rrij 


, V£ + 2<j<n-l, 






= m(f) 




(14) 



We can now repeat this kind of two-mode squeezing transformation between the {£ + 2)^^ 
mode and the -n}^ mode, and so on. Assume we are able to carry out this process only r 
times. We will have, at the end of it. 



m 



— i^mi, • • • , m^+r, m^_,_^_,_i, ■ ■ ■ , m„ 



mf^''^ = mf , V ^ + r + 1 < i < n - 1, 

m(f+'') = «n + Q-hi + Q-h2 H \- Q-hr, (15) 



so that = (5(0) - 2(e^+i + + ■ ■ ■ + e^+r)- Clearly, < < 2e^+,+i = 

2{me+r+i — "n^^^+i) (the last inequality encodes the fact that we could not carry out the 
Stage 2 operation one more time). 

Stage 3: Assume 5^^^'''^ > ( if 5^^+''^ = 0, we move directly to Stage 4). Carry out 
a two-mode canonical transformation between the (£ + r + l)**^ mode and the n^^ mode, 
taking the corresponding diagonal entries m^^i, ^^t'^ to m,+,+i = mf+;i + ^r+e+i and 
m(f+'^+-^) = m(f+^') + (5(^+^') — er+£+i respectively, leaving the other diagonals invariant, so that 
we have 

^(£+r+i) _ ^ Y£_|_,-_l_2<j<n — 1, 



n 

(^+r+l) 

j=£+r+2 j=£+r+2 



E mf'^-^^)= E (16) 



i.e., the situation in respect of the remaining n — {i+r + l) (or n — i — r if 5^^+'^) = 0) modes is 
precisely of the Schur-Horn type, suggesting that we deploy the beam splitter transformation 
n — I — r — 2 (or — / — r — 1 ) times. 

Stage 4 '■ Note that at the end of Stage 3 we have m^f"'"''"'"^) larger than m„ precisely 
by the sum of the amounts by which m^^'^^j-, for l<j<n — i — r — 2, are less than 
me+r+i+j- Therefore, for each value of j in this range, we effect a beam splitter transfor- 
mation connecting the {i + r + 1 + j)^^ mode to the n*^ mode, raising m^^'^^j to rrij and 
correspondingly pulling m^f"'"^'"''^^-') down by an equal amount. It is clear that at the end of 
these n — £ — r — 2 (or n — £ — r — l) redistributions, the diagonals will be precisely m. That 
is, m^""^-* = m. This completes proof of the theorem. 

We have taken maximal advantage of the simpler two- mode transformations Sg, S*^. The 
former was deployed r times in Stage 1 and n — £ — r — 2 (or n — £ — r — 1) times in Stage 4, 
and the latter £ times in Stage 2. The more general two-mode transformation was deployed 
(at the most) once in Stage 3. 



As illustration, and for comparison with Ref. we apply our procedure to the ex- 
ample noted after the statement of the theorem. The difference between I]J=i = 63 
and J2j=i i^j = 45 indicates the amount of squeezing that will have to be deployed at 
Stages 3 and 4. We have m(°) = k = (1,2,3,4,5,12,18); m« = (6,2,3,4,5,7,18); 
m(2) = (6,7,3,4,5,2,18); m^^) = (6,7,8,4,5,2,23); m^^) = (6,7,8,9,5,2,26); m^^) = 
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(6, 7, 8, 9, 10, 2, 21); and m(^) = (6, 7, 8, 9, 10, 11, 12) = m. The number of two-mode transfor- 
mations required at the four stages are 2, 1, 1, and 2 respectively. Note that m^'^^ dominates 
^(fc+i)^ for A; = 0,1, •••,5. 
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